The Mediator kinase module enhances polymerase activity to regulate transcriptional memory after heat stress in Arabidopsis.

Plants are often exposed to recurring adverse environmental conditions in the wild. Acclimation to high temperatures entails transcriptional responses, which prime plants to better withstand subsequent stress events. Heat stress (HS)-induced transcriptional memory results in more efficient re-induction of transcription upon recurrence of heat stress. Here, we identified CDK8 and MED12, two subunits of the kinase module of the transcription co-regulator complex, Mediator, as promoters of heat stress memory and associated histone modifications in Arabidopsis. CDK8 is recruited to heat-stress memory genes by HEAT SHOCK TRANSCRIPTION FACTOR A2 (HSFA2). Like HSFA2, CDK8 is largely dispensable for the initial gene induction upon HS, and its function in transcriptional memory is thus independent of primary gene activation. In addition to the promoter and transcriptional start region of target genes, CDK8 also binds their 3'-region, where it may promote elongation, termination, or rapid re-initiation of RNA polymerase II (Pol II) complexes during transcriptional memory bursts. Our work presents a complex role for the Mediator kinase module during transcriptional memory in multicellular eukaryotes, through interactions with transcription factors, chromatin modifications, and promotion of Pol II efficiency.

The gain-of-function mutation blf13 in the barley orthologue of the rice growth regulator NARROW LEAF1 is associated with increased leaf width.

Canopy architecture in cereals plays an important role in determining yield. Leaf width represents one key aspect of this canopy architecture. However, our understanding of leaf width control in cereals remains incomplete. Classical mutagenesis studies in barely identified multiple morphological mutants, including those with differing leaf widths. Of these, we characterized the broad leaf13 (blf13) mutant in detail. Mutant plants form wider leaves due to increased post-initiation growth and cell proliferation. The mutant phenotype perfectly co-segregated with a missense mutation in the HvHNT1 gene which affected a highly conserved region of the encoded protein, orthologous to the rice NARROW LEAF1 (NAL1) protein. Causality of this mutation for the blf13 phenotype is further supported by correlative transcriptomic analyses and protein-protein interaction studies showing that the mutant HvNHT1 protein interacts more strongly with a known interactor than wild-type HvHNT1. The mutant HvHNT1 protein also showed stronger homodimerization compared with wild-type HvHNT1, and homology modelling suggested an additional interaction site between HvHNT1 monomers due to the blf13 mutation. Thus, the blf13 mutation parallels known gain-of-function NAL1 alleles in rice that increase leaf width and grain yield, suggesting that the blf13 mutation may have a similar agronomic potential in barley.

Histone retention preserves epigenetic marks during heat stress-induced transcriptional memory in plants.

Plants often experience recurrent stressful events, for example, during heat waves. They can be primed by heat stress (HS) to improve the survival of more severe heat stress conditions. At certain genes, sustained expression is induced for several days beyond the initial heat stress. This transcriptional memory is associated with hyper-methylation of histone H3 lysine 4 (H3K4me3), but it is unclear how this is maintained for extended periods. Here, we determined histone turnover by measuring the chromatin association of HS-induced histone H3.3. Genome-wide histone turnover was not homogenous; in particular, H3.3 was retained longer at heat stress memory genes compared to HS-induced non-memory genes during the memory phase. While low nucleosome turnover retained H3K4 methylation, methylation loss did not affect turnover, suggesting that low nucleosome turnover sustains H3K4 methylation, but not vice versa. Together, our results unveil the modulation of histone turnover as a mechanism to retain environmentally mediated epigenetic modifications.

Testing the effect of individual scent compounds on pollinator attraction in nature using quasi-isogenic Capsella lines.

PREMISE: Floral scent, usually consisting of multiple compounds, is a complex trait, and its role in pollinator attraction has received increasing attention. However, disentangling the effect of individual floral scent compounds is difficult due to the complexity of isolating the effect of single compounds by traditional methods. METHODS: Using available quasi-isogenic lines (qILs) that were generated as part of the original mapping of the floral scent volatile-related loci CNL1 (benzaldehyde) and TPS2 (ß-ocimene) in Capsella, we generated four genotypes that should only differ in these two compounds. Plants of the four genotypes were introduced into a common garden outside the natural range of C. rubella or C. grandiflora, with individuals of a self-compatible C. grandiflora line as pollen donors, whose different genetic background facilitates the detection of outcrossing events. Visitors to flowers of all five genotypes were compared, and the seeds set during the common-garden period were collected for high-throughput amplicon-based sequencing to estimate their outcrossing rates. RESULTS: Benzaldehyde and ß-ocimene emissions were detected in the floral scent of corresponding genotypes. While some pollinator groups showed specific visitation preferences depending on scent compounds, the outcrossing rates in seeds did not vary among the four scent-manipulated genotypes. CONCLUSIONS: The scent-manipulated Capsella materials constructed using qILs provide a powerful system to study the ecological effects of individual floral scent compounds under largely natural environments. In Capsella, individual benzaldehyde and ß-ocimene emission may act as attractants for different types of pollinators.

MYB24 orchestrates terpene and flavonol metabolism as light responses to anthocyanin depletion in variegated grape berries.

Variegation is a rare type of mosaicism not fully studied in plants, especially fruits. We examined red and white sections of grape (Vitis vinifera cv. 'Béquignol') variegated berries and found that accumulation of products from branches of the phenylpropanoid and isoprenoid pathways showed an opposite tendency. Light-responsive flavonol and monoterpene levels increased in anthocyanin-depleted areas in correlation with increasing MYB24 expression. Cistrome analysis suggested that MYB24 binds to the promoters of 22 terpene synthase (TPS) genes, as well as 32 photosynthesis/light-related genes, including carotenoid pathway members, the flavonol regulator HY5 HOMOLOGUE (HYH), and other radiation response genes. Indeed, TPS35, TPS09, the carotenoid isomerase gene CRTISO2, and HYH were activated in the presence of MYB24 and MYC2. We suggest that MYB24 modulates ultraviolet and high-intensity visible light stress responses that include terpene and flavonol synthesis and potentially affects carotenoids. The MYB24 regulatory network is developmentally triggered after the onset of berry ripening, while the absence of anthocyanin sunscreens accelerates its activation, likely in a dose-dependent manner due to increased radiation exposure. Anthocyanins and flavonols in variegated berry skins act as effective sunscreens but for different wavelength ranges. The expression patterns of stress marker genes in red and white sections of 'Béquignol' berries strongly suggest that MYB24 promotes light stress amelioration but only partly succeeds during late ripening.

Genomic and epigenomic determinants of heat stress-induced transcriptional memory in Arabidopsis.

BACKGROUND: Transcriptional regulation is a key aspect of environmental stress responses. Heat stress induces transcriptional memory, i.e., sustained induction or enhanced re-induction of transcription, that allows plants to respond more efficiently to a recurrent HS. In light of more frequent temperature extremes due to climate change, improving heat tolerance in crop plants is an important breeding goal. However, not all heat stress-inducible genes show transcriptional memory, and it is unclear what distinguishes memory from non-memory genes. To address this issue and understand the genome and epigenome architecture of transcriptional memory after heat stress, we identify the global target genes of two key memory heat shock transcription factors, HSFA2 and HSFA3, using time course ChIP-seq. RESULTS: HSFA2 and HSFA3 show near identical binding patterns. In vitro and in vivo binding strength is highly correlated, indicating the importance of DNA sequence elements. In particular, genes with transcriptional memory are strongly enriched for a tripartite heat shock element, and are hallmarked by several features: low expression levels in the absence of heat stress, accessible chromatin environment, and heat stress-induced enrichment of H3K4 trimethylation. These results are confirmed by an orthogonal transcriptomic data set using both de novo clustering and an established definition of memory genes. CONCLUSIONS: Our findings provide an integrated view of HSF-dependent transcriptional memory and shed light on its sequence and chromatin determinants, enabling the prediction and engineering of genes with transcriptional memory behavior.

Poly(A) polymerase 1 contributes to competence acquisition of pollen tubes growing through the style in Arabidopsis thaliana.

Polyadenylation of mRNAs is critical for their export from the nucleus, stability, and efficient translation. The Arabidopsis thaliana genome encodes three isoforms of canonical nuclear poly(A) polymerase (PAPS) that redundantly polyadenylate the bulk of pre-mRNAs. However, previous studies have indicated that subsets of pre-mRNAs are preferentially polyadenylated by either PAPS1 or the other two isoforms. Such functional specialization raises the possibility of an additional level of gene-expression control in plants. Here we test this notion by studying the function of PAPS1 in pollen-tube growth and guidance. Pollen tubes growing through female tissue acquire the competence to find ovules efficiently and upregulate PAPS1 expression at the transcriptional, but not detectably at the protein level compared with in vitro grown pollen tubes. Using the temperature-sensitive paps1-1 allele we show that PAPS1 activity during pollen-tube growth is required for full acquisition of competence, resulting in inefficient fertilization by paps1-1 mutant pollen tubes. While these mutant pollen tubes grow almost at the wild-type rate, they are compromised in locating the micropyles of ovules. Previously identified competence-associated genes are less expressed in paps1-1 mutant than in wild-type pollen tubes. Estimating the poly(A) tail lengths of transcripts suggests that polyadenylation by PAPS1 is associated with reduced transcript abundance. Our results therefore suggest that PAPS1 plays a key role in the acquisition of competence and underline the importance of functional specialization between PAPS isoforms throughout different developmental stages.

Female self-incompatibility type in heterostylous Primula is determined by the brassinosteroid-inactivating cytochrome P450 CYP734A50.

Most flowering plants are hermaphrodites, with flowers having both male and female reproductive organs. One widespread adaptation to limit self-fertilization is self-incompatibility (SI), where self-pollen fails to fertilize ovules.1,2 In homomorphic SI, many morphologically indistinguishable mating types are found, although in heteromorphic SI, the two or three mating types are associated with different floral morphologies.3-6 In heterostylous Primula, a hemizygous supergene determines a short-styled S-morph and a long-styled L-morph, corresponding to two different mating types, and full seed set only results from intermorph crosses.7-9 Style length is controlled by the brassinosteroid (BR)-inactivating cytochrome P450 CYP734A50,10 yet it remains unclear what defines the male and female incompatibility types. Here, we show that CYP734A50 also determines the female incompatibility type. Inactivating CYP734A50 converts short S-morph styles into long styles with the same incompatibility behavior as L-morph styles, and this effect can be mimicked by exogenous BR treatment. In vitro responses of S- and L-morph pollen grains and pollen tubes to increasing BR levels could only partly explain their different in vivo behavior, suggesting both direct and indirect effects of the different BR levels in S- versus L-morph stigmas and styles in controlling pollen performance. This BR-mediated SI provides a novel mechanism for preventing self-fertilization. The joint control of morphology and SI by CYP734A50 has important implications for the evolutionary buildup of the heterostylous syndrome and provides a straightforward explanation for why essentially all of the derived self-compatible homostylous Primula species are long homostyles.11.

Phylogeography of a widely distributed plant species reveals cryptic genetic lineages with parallel phenotypic responses to warming and drought conditions.

To predict how widely distributed species will perform under future climate change, it is crucial to understand and reveal their underlying phylogenetics. However, detailed information about plant adaptation and its genetic basis and history remains scarce and especially widely distributed species receive little attention despite their putatively high adaptability. To examine the adaptation potential of a widely distributed species, we sampled the model plant Silene vulgaris across Europe. In a greenhouse experiment, we exposed the offspring of these populations to a climate change scenario for central Europe and revealed the population structure through whole-genome sequencing. Plants were grown under two temperatures (18°C and 21°C) and three precipitation regimes (65, 75, and 90 mm) to measure their response in biomass and fecundity-related traits. To reveal the population genetic structure, ddRAD sequencing was employed for a whole-genome approach. We found three major genetic clusters in S. vulgaris from Europe: one cluster comprising Southern European populations, one cluster of Western European populations, and another cluster containing central European populations. Population genetic diversity decreased with increasing latitude, and a Mantel test revealed significant correlations between F ST and geographic distances as well as between genetic and environmental distances. Our trait analysis showed that the genetic clusters significantly differed in biomass-related traits and in the days to flowering. However, half of the traits showed parallel response patterns to the experimental climate change scenario. Due to the differentiated but parallel response patterns, we assume that phenotypic plasticity plays an important role for the adaptation of the widely distributed species S. vulgaris and its intraspecific genetic lineages.

Heteromeric HSFA2/HSFA3 complexes drive transcriptional memory after heat stress in Arabidopsis.

Adaptive plasticity in stress responses is a key element of plant survival strategies. For instance, moderate heat stress (HS) primes a plant to acquire thermotolerance, which allows subsequent survival of more severe HS conditions. Acquired thermotolerance is actively maintained over several days (HS memory) and involves the sustained induction of memory-related genes. Here we show that FORGETTER3/ HEAT SHOCK TRANSCRIPTION FACTOR A3 (FGT3/HSFA3) is specifically required for physiological HS memory and maintaining high memory-gene expression during the days following a HS exposure. HSFA3 mediates HS memory by direct transcriptional activation of memory-related genes after return to normal growth temperatures. HSFA3 binds HSFA2, and in vivo both proteins form heteromeric complexes with additional HSFs. Our results indicate that only complexes containing both HSFA2 and HSFA3 efficiently promote transcriptional memory by positively influencing histone H3 lysine 4 (H3K4) hyper-methylation. In summary, our work defines the major HSF complex controlling transcriptional memory and elucidates the in vivo dynamics of HSF complexes during somatic stress memory.

Fairy circles in Namibia are assembled from genetically distinct grasses.

Fairy circles are striking regularly sized and spaced, bare circles surrounded by Stipagrostis grasses that occur over thousands of square kilometres in Namibia. The mechanisms explaining their origin, shape, persistence and regularity remain controversial. One hypothesis for the formation of vegetation rings is based on the centrifugal expansion of a single individual grass plant, via clonal growth and die-back in the centre. Clonality could explain FC origin, shape and long-term persistence as well as their regularity, if one clone competes with adjacent clones. Here, we show that for virtually all tested fairy circles the periphery is not exclusively made up of genetically identical grasses, but these peripheral grasses belong to more than one unrelated genet. These results do not support a clonal explanation for fairy circles. Lack of clonality implies that a biological reason for their origin, shape and regularity must emerge from competition between near neighbor individuals within each fairy circle. Such lack of clonality also suggests a mismatch between longevity of fairy circles versus their constituent plants. Furthermore, our findings of lack of clonality have implications for some models of spatial patterning of fairy circles that are based on self-organization.

Supergene evolution via stepwise duplications and neofunctionalization of a floral-organ identity gene.

Heterostyly represents a fascinating adaptation to promote outbreeding in plants that evolved multiple times independently. While l-morph individuals form flowers with long styles, short anthers, and small pollen grains, S-morph individuals have flowers with short styles, long anthers, and large pollen grains. The difference between the morphs is controlled by an S-locus "supergene" consisting of several distinct genes that determine different traits of the syndrome and are held together, because recombination between them is suppressed. In Primula, the S locus is a roughly 300-kb hemizygous region containing five predicted genes. However, with one exception, their roles remain unclear, as does the evolutionary buildup of the S locus. Here we demonstrate that the MADS-box GLOBOSA2 (GLO2) gene at the S locus determines anther position. In Primula forbesii S-morph plants, GLO2 promotes growth by cell expansion in the fused tube of petals and stamen filaments beneath the anther insertion point; by contrast, neither pollen size nor male incompatibility is affected by GLO2 activity. The paralogue GLO1, from which GLO2 arose by duplication, has maintained the ancestral B-class function in specifying petal and stamen identity, indicating that GLO2 underwent neofunctionalization, likely at the level of the encoded protein. Genetic mapping and phylogenetic analysis indicate that the duplications giving rise to the style-length-determining gene CYP734A50 and to GLO2 occurred sequentially, with the CYP734A50 duplication likely the first. Together these results provide the most detailed insight into the assembly of a plant supergene yet and have important implications for the evolution of heterostyly.

Mapping-by-Sequencing via MutMap Identifies a Mutation in ZmCLE7 Underlying Fasciation in a Newly Developed EMS Mutant Population in an Elite Tropical Maize Inbred.

Induced point mutations are important genetic resources for their ability to create hypo- and hypermorphic alleles that are useful for understanding gene functions and breeding. However, such mutant populations have only been developed for a few temperate maize varieties, mainly B73 and W22, yet no tropical maize inbred lines have been mutagenized and made available to the public to date. We developed a novel Ethyl Methanesulfonate (EMS) induced mutation resource in maize comprising 2050 independent M2 mutant families in the elite tropical maize inbred ML10. By phenotypic screening, we showed that this population is of comparable quality with other mutagenized populations in maize. To illustrate the usefulness of this population for gene discovery, we performed rapid mapping-by-sequencing to clone a fasciated-ear mutant and identify a causal promoter deletion in ZmCLE7 (CLE7). Our mapping procedure does not require crossing to an unrelated parent, thus is suitable for mapping subtle traits and ones affected by heterosis. This first EMS population in tropical maize is expected to be very useful for the maize research community. Also, the EMS mutagenesis and rapid mapping-by-sequencing pipeline described here illustrate the power of performing forward genetics in diverse maize germplasms of choice, which can lead to novel gene discovery due to divergent genetic backgrounds.

A Similar Genetic Architecture Underlies the Convergent Evolution of the Selfing Syndrome in Capsella.

Whether, and to what extent, phenotypic evolution follows predictable genetic paths remains an important question in evolutionary biology. Convergent evolution of similar characters provides a unique opportunity to address this question. The transition to selfing and the associated changes in flower morphology are among the most prominent examples of repeated evolution in plants. In this study, we take advantage of the independent transitions to self-fertilization in the genus Capsella to compare the similarities between parallel modifications of floral traits and test for genetic and developmental constraints imposed on flower evolution in the context of the selfing syndrome. Capsella rubella and Capsella orientalis emerged independently but evolved almost identical flower characters. Not only is the evolutionary outcome identical but the same developmental strategies underlie the convergent reduction of flower size. This has been associated with convergent evolution of gene expression changes. The transcriptomic changes common to both selfing lineages are enriched in genes with low network connectivity and with organ-specific expression patterns. Comparative genetic mapping also suggests that, at least in the case of petal size evolution, these similarities have a similar genetic basis. Based on these results, we hypothesize that the limited availability of low-pleiotropy paths predetermines closely related species to similar evolutionary outcomes.

Proteomic and metabolomic profiling underlines the stage- and time-dependent effects of high temperature on grape berry metabolism.

Climate change scenarios predict an increase in mean air temperatures and in the frequency, intensity, and length of extreme temperature events in many wine-growing regions worldwide. Because elevated temperature has detrimental effects on berry growth and composition, it threatens the economic and environmental sustainability of wine production. Using Cabernet Sauvignon fruit-bearing cuttings, we investigated the effects of high temperature (HT) on grapevine berries through a label-free shotgun proteomic analysis coupled to a complementary metabolomic study. Among the 2,279 proteins identified, 592 differentially abundant proteins were found in berries exposed to HT. The gene ontology categories "stress," "protein," "secondary metabolism," and "cell wall" were predominantly altered under HT. High temperatures strongly impaired carbohydrate and energy metabolism, and the effects depended on the stage of development and duration of treatment. Transcript amounts correlated poorly with protein expression levels in HT berries, highlighting the value of proteomic studies in the context of heat stress. Furthermore, this work reveals that HT alters key proteins driving berry development and ripening. Finally, we provide a list of differentially abundant proteins that can be considered as potential markers for developing or selecting grape varieties that are better adapted to warmer climates or extreme heat waves.

Retracing the molecular basis and evolutionary history of the loss of benzaldehyde emission in the genus Capsella.

The transition from pollinator-mediated outbreeding to selfing has occurred many times in angiosperms. This is generally accompanied by a reduction in traits attracting pollinators, including reduced emission of floral scent. In Capsella, emission of benzaldehyde as a main component of floral scent has been lost in selfing C. rubella by mutation of cinnamate-CoA ligase CNL1. However, the biochemical basis and evolutionary history of this loss remain unknown, as does the reason for the absence of benzaldehyde emission in the independently derived selfer Capsella orientalis. We used plant transformation, in vitro enzyme assays, population genetics and quantitative genetics to address these questions. CNL1 has been inactivated twice independently by point mutations in C. rubella, causing a loss of enzymatic activity. Both inactive haplotypes are found within and outside of Greece, the centre of origin of C. rubella, indicating that they arose before its geographical spread. By contrast, the loss of benzaldehyde emission in C. orientalis is not due to an inactivating mutation in CNL1. CNL1 represents a hotspot for mutations that eliminate benzaldehyde emission, potentially reflecting the limited pleiotropy and large effect of its inactivation. Nevertheless, even closely related species have followed different evolutionary routes in reducing floral scent.

A high-throughput amplicon-based method for estimating outcrossing rates.

BACKGROUND: The outcrossing rate is a key determinant of the population-genetic structure of species and their long-term evolutionary trajectories. However, determining the outcrossing rate using current methods based on PCR-genotyping individual offspring of focal plants for multiple polymorphic markers is laborious and time-consuming. RESULTS: We have developed an amplicon-based, high-throughput enabled method for estimating the outcrossing rate and have applied this to an example of scented versus non-scented Capsella (Shepherd's Purse) genotypes. Our results show that the method is able to robustly capture differences in outcrossing rates. They also highlight potential biases in the estimates resulting from differential haplotype sharing of the focal plants with the pollen-donor population at individual amplicons. CONCLUSIONS: This novel method for estimating outcrossing rates will allow determining this key population-genetic parameter with high-throughput across many genotypes in a population, enabling studies into the genetic determinants of successful pollinator attraction and outcrossing.

The poly(A) polymerase PAPS1 interacts with the RNA-directed DNA-methylation pathway in sporophyte and pollen development.

RNA-based processes play key roles in the regulation of eukaryotic gene expression. This includes both the processing of pre-mRNAs into mature mRNAs ready for translation and RNA-based silencing processes, such as RNA-directed DNA methylation (RdDM). Polyadenylation of pre-mRNAs is one important step in their processing and is carried out by three functionally specialized canonical nuclear poly(A) polymerases in Arabidopsis thaliana. Null mutations in one of these, termed PAPS1, result in a male gametophytic defect. Using a fluorescence-labelling strategy, we have characterized this defect in more detail using RNA and small-RNA sequencing. In addition to global defects in the expression of pollen-differentiation genes, paps1 null-mutant pollen shows a strong overaccumulation of transposable element (TE) transcripts, yet a depletion of 21- and particularly 24-nucleotide-long short interfering RNAs (siRNAs) and microRNAs (miRNAs) targeting the corresponding TEs. Double-mutant analyses support a specific functional interaction between PAPS1 and components of the RdDM pathway, as evident from strong synergistic phenotypes in mutant combinations involving paps1, but not paps2 paps4, mutations. In particular, the double-mutant of paps1 and rna-dependent rna polymerase 6 (rdr6) shows a synergistic developmental phenotype disrupting the formation of the transmitting tract in the female gynoecium. Thus, our findings in A. thaliana uncover a potentially general link between canonical poly(A) polymerases as components of mRNA processing and RdDM, reflecting an analogous interaction in fission yeast.

Successive duplication-divergence mechanisms at the RCO locus contributed to leaf shape diversity in the Brassicaceae.

Gene duplication is a major driver for the increase of biological complexity. The divergence of newly duplicated paralogs may allow novel functions to evolve, while maintaining the ancestral one. Alternatively, partitioning the ancestral function among paralogs may allow parts of that role to follow independent evolutionary trajectories. We studied the REDUCED COMPLEXITY (RCO) locus, which contains three paralogs that have evolved through two independent events of gene duplication, and which underlies repeated events of leaf shape evolution within the Brassicaceae. In particular, we took advantage of the presence of three potentially functional paralogs in Capsella to investigate the extent of functional divergence among them. We demonstrate that the RCO copies control growth in different areas of the leaf. Consequently, the copies that are retained active in the different Brassicaceae lineages contribute to define the leaf dissection pattern. Our results further illustrate how successive gene duplication events and subsequent functional divergence can increase trait evolvability by providing independent evolutionary trajectories to specialized functions that have an additive effect on a given trait.

Variation in Splicing Efficiency Underlies Morphological Evolution in Capsella.

Understanding the molecular basis of morphological change remains a central challenge in evolutionary-developmental biology. The transition from outbreeding to selfing is often associated with a dramatic reduction in reproductive structures and functions, such as the loss of attractive pheromones in hermaphroditic Caenorhabditis elegans and a reduced flower size in plants. Here, we demonstrate that variation in the level of the brassinosteroid-biosynthesis enzyme CYP724A1 contributes to the reduced flower size of selfing Capsella rubella compared with its outbreeding ancestor Capsella grandiflora. The primary transcript of the C. rubella allele is spliced more efficiently than that of C. grandiflora, resulting in higher brassinosteroid levels. These restrict organ growth by limiting cell proliferation. More efficient splicing of the C. rubella allele results from two de novo mutations in the selfing lineage. Thus, our results highlight the potentially widespread importance of differential splicing efficiency and higher-than-optimal hormone levels in generating phenotypic variation.